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Abstract 

Liolaemus audituvelatus (Nunez & Yanez 1983) and L. manueli (Nunez, Navarro, Garin, Pincheira-Donoso & Meriggio 
2003) are endemic species of the Atacama Desert of northern Chile that belong to the montanus group. Both species are 
considered cryptic from each other and can only be distinguished by their distribution ranges and karyotypes. Originally, 
there was a wide separation zone between their known distribution ranges, but later collections reduced the gap from 430 
km to only 150 km. In this study, we review the geographic information about both species and report new localities within 
the distribution gap, where species identification becomes difficult. We performed a molecular phylogenetic analysis and 
applied several species delimitation methods to reassess the taxonomic status of both nominal species and new interme¬ 
diate populations. Our analyses support the placement of L. manueli in the synonymy of L. audituvelatus. We discuss the 
biogeographic and conservation implications of this new synonymy. 
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Introduction 

Liolaemus is a lizard genus widely distributed throughout South America, ranging from central Peru and 
southeastern Brazil to Tierra del Fuego in Argentina and Chile (Pincheira-Donoso el al. 2008). The more than 250 
species of the genus (Abdala & Quinteros 2014; Uetz & Hosek 2016) occupy the most diverse environments with 
altitudes ranging from sea level to approximately 5000 m in the Andean highlands (Donoso-Barros 1966; Cei 
1993; Pincheira-Donoso et al. 2008; Aparicio & Ocampo 2010). 

Two major clades have been recognized within Liolaemus'. The Chilean group proposed as the subgenus 
Liolaemus, and the Argentinian group, proposed as the subgenus Eulaemus (Laurent 1983; 1985; Etheridge 1995; 
Schulte et al. 2000). A subset of species belonging to the Argentinian group have been identified as the montanus 
group (Etheridge 1995), where the knowledge about many species is still poor and where there have been some 
taxonomic problems (Nunez & Jaksic 1992; Pincheira-Donoso & Nunez 2002; Valladares et al. 2002; Pincheira- 
Donoso & Nunez 2005; Troncoso-Palacios 2014; Ruiz de Gamboa & Ortiz 2016; Aguilar et al. 2016). This is 
probably because these species inhabit the desert or highlands where is difficult to study them. 

Within this group, two species of similar morphology were described in northern Chile: Liolaemus 
audituvelatus (Nunez & Yanez 1983) and L. manueli (Nunez, Navarro, Garin, Pincheira-Donoso & Meriggio 
2003). Originally, L. audituvelatus was known only in the type locality Llano de Vilama, 10 km east of San Pedro 
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de Atacama (Nunez & Yanez 1983). Years later, Nunez et al. (2003) mentioned that the distribution of L. 
audituvelatus was restricted to the Salar de Atacama (50 km south of San Pedro de Atacama) in Antofagasta 
Region. Liolaemus manueli was known only from its type locality: Diego de Almagro in the Antofagasta Region 
(Nunez et al. 2003). Between their type localities, there existed a wide separation zone (430 km in straight line), 
where it was assumed that there were no populations connecting both species. Differentiated morphologically from 
each other only by imprecise subtleties, they were considered cryptic species, distinguishable only by their 
distribution and karyotypes (Nunez et al. 2003). 

Later, Nunez et al. (2012) reported new localities for both species: L. manueli was recorded 124 km southwest 
of its type locality and L. audituvelatus was recorded 280 km southwest of its type locality. This last record reduced 
the gap separating the known localities of the two species, from 430 km to only 150 km. Due to the lack of 
diagnostic morphological characters between these species, the correct assignment of individuals to species is 
complicated in intermediate records. In this study, we update the geographic data of both species and evaluate 
whether L. audituvelatus and L. manueli correspond to different phylogenetic taxa by applying a series of species 
delimitation methods. 


Materials and methods 

For this study, we compiled all the geographic records of the literature and the database of the Museo Nacional de 
Historia Natural de Chile (MNHN) herpetological collection (GBIF, 2016), to which we added new undescribed 
localities (Figure 1, Appendix SI) assigned to Liolaemus audituvelatus (Figure 2a) or L. manueli by geographic 
criteria. We also collected specimens near the type locality of L. audituvelatus (near Caspana), and other 
intermediate localities between the known distributions of L. manueli (Figures 2d, e and f) and L. audituvelatus 
(Figure 1) (Figure 2b and c). Additionally, a total of 16 individuals belonging to seven species of the montanus 
group were included in the phylogenetic and species delimitation analyses (Appendix S2). 

The collected individuals were sacrificed by injection of sodium thiopental, applied in the pericardial cavity. 
For the molecular study, we cut a small piece of the tail or liver, which were immediately stored in absolute ethanol 
(>99.5%). 

Genetic data and analyses. Genomic DNA was extracted from the tail or liver using a Promega Wizard SV 
Genomic DNA Purification System Kit, following the manufacturer’s protocol. A fragment of the mitochondrial 
gene cytochrome b (Kocher et al. 1989) was amplified using the PCR protocol of Vidal et al. (2012) and sequenced 
with an Applied Biosystems 3730XL sequencer (Macrogen, Korea). Sequences were edited with CodonCode 
Aligner software (CodonCode Corporation, Dedham, MA.) and translated to amino acids to identify open reading 
frames. 

As an ingroup, we used specimens that could be assigned to L. audituvelatus by distribution range, including 
one found near the type locality (Figure 3), one specimen from the type locality of/., manueli and another one from 
a known locality, specimens of intermediate populations between the type localities of these two species, and 
representatives of other species of the Liolaemus ( Eulaemus ) montanus group (Schulte et al. 2000) from Chile and 
another for which cytb sequences are available in GenBank (Appendix S2). As an outgroup, we used L. ornatus 
Koslowsky, belonging to the Liolaemus (Eulaemus) darwinii group (Aguilar et al. 2016). All new sequences were 
deposited in GenBank (Appendix S2). 

Sequences were aligned with MUSCLE (Edgar 2004). We evaluated whether the sequences were saturated and 
thus useful for the phylogenetic analysis, using the Xia test (Xia et al. 2003; Xia & Lemey 2009) implemented in 
DAMBE v6.4.29 (Xia 2013). This is an entropy-based index that estimates a substitution saturation index (Iss) and 
compares it to a critical substitution saturation index (Iss.c) via a randomization process with 95% confidence 
intervals (Xia 2000; Xia & Lemey 2009). Phylogenetic relationships were inferred using a Markov Chain Monte 
Carlo (MCMC) method within a Bayesian framework (BMCMC) to estimate the posterior probability of 
phylogenetic trees in BayesPhylogenies vl.l software (http://www.evolution.rdg.ac.uk/ BayesPhy.html). We ran 
BMCMC analyses consisting of 10,000,000 iterations of phylogenetic trees, sampled every 1000th to assure that 
successive samples were independent. From this sample of trees, we removed the first thousand to avoid including 
samples before the convergence of the Markov Chains. 
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FIGURE 1. Distribution maps of Liolaemits audituvelatus and L. manueli based on historical and new records. Type localities 
of both species are indicated by stars. Historical records of both species gathered from Nunez et al. (1998, 2003, 2012) and 
GBIF (2016) also are indicated by black circles and squares. New localities of L. audituvelatus are indicated by red circles and 
red squares for L. manueli. Localities from which DNA samples were obtained for the analyses of this study are indicated with 
green triangles. The geographic references are in Appendix S2. 

Species delimitation. We applied three species delimitation methods to evaluate whether L. audituvelatus and 
L. manueli are phylogenetically distinct. We used Automatic Barcode Gap Discovery (ABGD) (Puillandre et al. 
2012), Poisson Tree Process (PTP) (Zhang et al. 2013) and multi-rate Poisson Tree Process (mPTP) (Kapli et al. 
2017), which are based on different aspects of molecular evolution. The Automatic Barcode Gap Discovery 
(ABGD) method was used to divide data into candidate species based on a statistically inferred barcode gap 
(Puillandre et al. 2012), employing three distance options (Simple, Jukes-Cantor and Kimura 2-parameter). We 
used a relative gap width (X) of 1.0 and tested a range between 0.001 and 0.25 of prior intraspecific divergence 
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values (Puillandre el al. 2012). The corresponding analysis was conducted on the web server for ABGD (http:// 
wwwabi.snv.jussieu.fr/public/abgd/abgdweb.html). The PTP method models speciation rate by directly using the 
number of substitutions and applying two independent Poisson process classes, one which describes speciation and 
another which describes the coalescent intraspecific processes. Thus, the adjustment of both processes delimits the 
species in a given topology (Zhang et al. 2013). A new method based on the PTP algorithm, multi-rate PTP 
(mPTP), improves the estimate considering different rates of speciation-coalescence (Kapli et al. 2016). Although 
the speciation rate can be assumed as constant between closely related species, intraspecific coalescence rate and 
genetic diversity can differ significantly even among sister species, then mPTP allows to consider different rates of 
branching events within each delimited species (Kapli et al. 2016). To perform PTP and mPTP we used the web 
server (http://mptp.h-its.org/) and the Bayesian consensus tree. For all these analyses, the sequence of the outgroup 
was excluded. 

Moreover, we complemented the species delimitation analyses by computing pairwise distances among all 
specimens included in the phylogenetic analysis with BioEdit v.7.0.9 (Flail 1999). 



FIGURE 2. Individuals from a) near Caspana (MUAP-114), near type locality for Liolaemus audituvelatus; b) Diego de 
Almagro, type locality of L. manueli; c) and d) Altos Quebrada Agua Colorada; e) and f) Barranquilla (near Caseron). fj is a 
gravid female. Note that b) male of L. audituvelatus from known distribution and d) male of L. manueli from type locality, are 
identical and were found 300 km away (in straight line); and, a) and e) are also similar and were found 600 km away (in straight 
line). 


Results 

Phylogenetic Analysis. The Bayesian analysis included cytb sequences from a total of 29 specimens assigned to 
15 nominal species of the montanus group, seven from Argentina and eight from Chile. The consensus tree with 
Bayesian posterior probabilities, based on 705 nucleotide sites, is shown in Figure 4. The analysis recovered L. 
audituvelatus and L. manueli as a monophyletic group with high support (BPP = 1), but nevertheless, the two 
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specimens of L. audituvelatus from near Caspana are paraphyletic with respect to all samples assigned to L. 
manueli. 



FIGURE 3. Specimens of Liolaemus audituvelatus collected near Caspana, 65 km north of San Pedro de Atacama, 
Antofagasta Region (MUAP-113 (LNC135) and MUAP-114 (LNC136)). 
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FIGURE 4. Bayesian consensus tree showing the phylogenetic relationships among species of the montanus group and new 
populations related to L. audituvelatus and /,. manueli. Posterior probability is indicated at each node. 
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Species delimitation analyses. Each species delimitation method recognized a different number of candidate 
species (Figure 5). ABGD recognized 11 candidate species within each of the three distances used. This method 
suggests that L. audituvelatus and L. manueli correspond to a single species, as was also the case for the clade 
including L. eleodori Cei, Etheridge & Videla, L. vallecurensis Pereyra, L. rosenmanni Nunez & Navarro and L. 
nigriceps (Philippi). 

The mPTP model recognized 10 candidate species and also suggested that L. audituvelatus and L. manueli 
correspond to a single species. Likewise, the clade including L. eleodori , L. vallecurensis, L. rosenmanni and L. 
nigriceps was also found to correspond to a single species as was the clade including L. andinus Koslowsky and L. 
famatinae Cei. 

The PTP model recognized 14 candidate species. This method also suggests that L. audituvelatus and L. 
manueli are conspecifics, but one individual of L. audituvelatus (MUAP-114 (LNC136) was recognized as a 
different species. Also, L. eleodori with L. vallecurensis likely correspond to the same species. 

While each of the three species delimitation methods found Liolaemus audituvelatus and L. manueli to belong 
to a single species, one of the two specimens from near Caspana (MUAP-114 (LNC136); Figure 3) was assigned to 
L. audituvelatus was recognized as a separate candidate species by the mPTP analysis. All methods placed the 
other Caspana specimen (MUAP-113 (LNC135)) in the same species as all samples assigned to L. manueli (Figure 
5). It is necessary to emphasize that both specimens of L. audituvelatus collected near Caspana are 
indistinguishable morphologically (Figure 3) and that the genetic distance between specimen MUAP-113 
(LNC135) and the specimen from the type locality of L. manueli (Diego of Almagro) is extremely low (0.0043, 
Appendix S3), even though the localities are 480 km apart. 

Genetic distances. Considering the current nominal taxonomy, intraspecies genetic distances fluctuated 
between 0 and 0.0172 (Appendix S3). The highest value corresponded to the two individuals of Liolaemus 
audituvelatus that were collected near Caspana. The second highest value (0.0143) corresponds to the two 
individuals of L. torresi (Nunez, Navarro, Garin, Pincheira-Donoso & Meriggio), which were collected 60 km 
apart. Lower intraspecies genetic distances were obtained among the samples of L. manueli, with a maximum value 
of 0.0082 between two individuals collected 120 km apart. 

Interspecies genetic distances with the sampled species of the montanus group fluctuated between 0 and 
0.1407 (Appendix S3). The highest distance was observed between L. signifer (Dumeril & Bibron) and L. 
poconchilensis Valladares, whereas L. eleodori and L. vallecurensis exhibited the lowest distance according to 
sequences generated in other studies (Morando et al. 2003; Olave el al. 2014). Liolaemus audituvelatus and L. 
manueli have the second lowest distance (0.0043) between individuals which were collected at 480 km away. 
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FIGURE 5. Results of the species delimitation analyses. On the left, the Bayesian consensus topology used for PTP and mPTP 
analyses is shown. On the right, boxes represent the species delimitations obtained with the three methods, in comparison with 
the current taxonomy (column nominal taxa). 
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Discussion 


Although the delimitation methods we employed differ in the number of candidate species recognized, they are 
congruent enough to suggest that L. audituvelatus and L. manueli belong to a single species. The discrepancies 
among analyses are few and are found in non-focal species. For example, the clade ( L. eleodori + L. vallecurensis ) 
+ ( L. rosenmanni + L. nigriceps) was recognized as one species in the ABGD and mPTP analyses. Species of this 
clade inhabit Argentina and Chile and exhibit enough morphological differences (Philippi 1860, Cei et al. 1985, 
Nunez & Navarro 1992, Pereyra 1992) to have been considered valid taxa (Pincheira-Donoso & Nunez 2005; 
Pincheira-Donoso et al. 2008; Lobo et al. 2010; Abdala & Quinteros 2014), but they are separated by low genetic 
distances (< 2%) according to all available data (Morando et al. 2003; Olave et al. 2014; this study). Nevertheless, 
there exists a possibility that some of these species are not valid because L. rosenmanni was initially reported as Z. 
eleodori in Chile (Nunez & Torres-Mura 1992; Nunez & Jaksic 1992; Nunez & Navarro 1992). Moreover, it would 
be useful to review the identification of the individuals recognized as L. vallecurensis by Morando et al. (2003) and 
L. eleodori by Olave et al. (2014) because they display low genetic differences in cytochrome b (Appendix S3). 
Consequently, our delimitation analyses recovered them as conspecifics. These results suggest that more in-depth 
studies are needed, specifically studies including more specimens that also expand the geographical distribution of 
the samples in order to morphologically and phylogenetically delimit the valid taxa among these species. 

Although mPTP is the most sophisticated methodology used here, in our case it underestimates the number of 
species when current taxonomy is taken into account. PTP was the method that most coincided with the currently 
accepted taxonomy. 

The mPTP analysis recognized both specimens collected near Caspana to be belong to different species, this 
suggests that there are two cryptic species at Caspana, with one restricted to that locality and the other with a wide 
distribution towards the south. This would imply the sympatry of two sister species in one locality, without evident 
morphological differences, which poses a more complex and improbable evolutionary scenario. Moreover, it is 
noteworthy that the other two analyses (ABGD and PTP) support that they are a single species together with all 
samples assigned to L. manueli. We infer that there would a higher intraspecific genetic divergence in the 
population of Caspana which now defines the northern end of the distribution of L. audituvelatus. 

So far, only the karyotype allows for the differentiation of L. audituvelatus from L. manueli. Nevertheless, 
chromosomes may vary in form, size and even in number, within and between populations of Liolaemus (Lamborot 
et al. 1979; 1981; 1993; 2006; Lamborot & Alvarez-Sarret 1989; 1993; Lamborot 1991; 1998; 2001; Iturra et al. 
1994; Lamborot & Vasquez 1998); even the descriptions of karyotypes of the same species may differ when 
comparing different sources (e.g. Navarro Baron 1991; Aiassa et al. 1999; Aiassa et al. 2005, Aiassa & Gorla 
2010). Due to this variability in the karyotype as a character, this would not be reliable enough to define species in 
Liolaemus and therefore is not enough to differentiate to L. manueli of L. audituvelatus. 

Liolaemus audituvelatus was described in the genus Ctenoblepharys (Nunez & Yanez 1983) and L. manueli in 
Phrynosaura (Nunez et al. 2003). Later, C. audituvelatus was incorporated into Phrynosaura (Laurent 1984; 
Nunez & Yanez 1984). Although Phrynosaura was synonymized with Liolaemus more than twenty years ago 
(Frost & Etheridge 1989; Etheridge 1995), several authors have continued considering it a valid genus (Nunez & 
Jaksic 1992; Nunez et al. 1998; 2003; Pincheira-Donoso & Nunez 2005; Vidal et al. 2009) or a group of species 
within Liolaemus ( L. reichei group in Pincheira-Donoso et al. 2008). Our phylogenetic hypothesis includes all 
species described or included in Phrynosaura (Z. audituvelatus, L. manueli, L. poconchilensis, L. stolzmanni and Z. 
torresi) except Z. erroneus, which was described by a single specimen, now lost, and whose type locality is not 
known with exactitude. For this reason, no later records exist. Our results support and confirms that Phrynosaura 
or the Z. reichei group are not monophyletic (Valladares et al. 2002; Lobo et al. 2010), because that the species 
formerly recognized by some authors as belonging to Phrynosaura are placed in separate clades by the present 
phylogenetic hypothesis. 

Based on their great morphological similarity and their close phylogenetic affinity, we consider Z. 
audituvelatus and Z. manueli to be conspecific and here recognize Z. manueli to be a junior synonym of Z. 
audituvelatus. Our placement of Liolaemus manueli in the synonymy of Z. audituvelatus implies a considerable 
range extension of the senior species, which affects its conservation status. Liolaemus audituvelatus is categorized 
as Vulnerable Blab(iii); D2 by the IUCN (Nunez et al. 2017) and is considered Vulnerable Blab(iii)+2ab(iii) by 
Chilean Law (D.S. N° 16/2016 MMA), while Z. manueli is considered Endangered Blab(iii)+2ab(iii) in Chile 
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(D.S. N° 16/2016 MMA) and Endangered Blab(iii) by the IUCN (Valladares el al. 2016). It is worth noting that L. 
audituvelatus is a desert species, living in places of very scant to no vegetation and populations of low densities. 
Thus, it is hard to observe in the field. Moreover, the places where it lives are under pressure from intensive mining 
activities and vehicular traffic due to recreational activities like raids. Based on the considered range expansion of 
this taxon, we recommend that L. audituvelatus , in its new definition, remains Vulnerable. 
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MUAP-104 MH184782 L.manueli Llano de Barranquilla. Atacama -27.3803 -70.6554 YM1 

Region 

MUAP-118 MH184781 L.manueli Diego de Almagro. Atacama Region -26.4733 -69.9664 YM2 

MUAP-110 MH184779 L.manueli Diego de Alamgro. km 65 ruta C-13. -26.35357 -69.92448 LNC86 











APPENDIX S2. (Continued) 

VOUCHER GENBANK SPECIES LOCALITY LATITUDE LONGITUDE OBSERVATION 

_ ACCESSION _ 

SSUC-Re 760 MH184783 L. manueli Llano La Polvora. Antofagasta Region -25.81741 -69.527632 ER2 

ER1 MH184780 L. manueli 17 km S of Las Breas. Ruta A-5 near -25.71352 -70.354936 Released on the same 
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MUAP-116 MH184797 L.torresi near Maria Elena. Antofagasta Region LNC146 

AY173546 L. vallecurensis Morando e? al. 2003 






APPENDIX S3. Pairwise genetic distances among the specimens used in the phylogenetic/delimitation analyses. 
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APPENDIX S3. (Continued) 

N° SEQUENCE _ 
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